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ABSTRACT 

We employ the first fully three-dimensional simulation to study the role of magnetic 
fields and ion-neutral friction in regulating gravitationally-driven fragmentation of 
molecular clouds. The cores in an initially subcritical cloud develop gradually over an 
ambipolar diffusion time while the cores in an initially supercritical cloud develop in 
a dynamical time. The infall speeds on to cores are subsonic in the case of an initially 
subcritical cloud, while an extended (> 0.1 pc) region of supersonic infall exists in 
the case of an initially supercritical cloud. These results are consistent with previous 
two-dimensional simulations. We also found that a snapshot of the relation between 
density (p) and the strength of the magnetic field (B) at different spatial points of the 
cloud coincides with the evolutionary track of an individual core. When the density 
becomes large, both relations tend to B oc p 05 . 

Key words: instabilities - ISM: clouds - ISM: magnetic fields - ISM: molecules - 
MHD - stars: formation 
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1 INTRODUCTION 

Magnetic fields in molecular clouds play an important role 
in the early stages of star formation. They may regulate 
the cloud fragmentation process, moderate the infall mo- 
tions on to density peaks, control angular momentum evo- 
lution through magnetic braking, launch jets from the near- 
protostellar environment, and possibly determine a finite 
mass reservoir for star formation by limiting accretion from a 
magnetically-dominated envelope. The prevailing macrotur- 
bulence in molecular cloud envelopes also likely represents 
magnetohydrodynamic motions. This paper concerns itself 
with the first two issues above: we employ the first fully 
three-dimensional simulation to study the role of magnetic 
fields and ion-neutral friction in regulating gravitationally- 
driven fragmentation and infall. 

The relative strengths of gravity and the mag- 
netic field can be quantified through the mass- 
to-flux ratio M/$. There exists a critica l mass- 
to-flux ratio CM/$ Uit. jMestel fc Spitzerl Il956 
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M/<& > (A//$) C rit, a pressure-bounded cloud is super- 
critical and is prone to indefinite collapse if the external 
pressure exceeds some value. This behaviour is analogous to 
that of the nonmagnetic Bonnor-Ebert sphere. Conversely, 
if M/$ < (M/$) c rit, a cloud is subcritical and cannot 
collapse even in the limit of infinite external pressure, as 
long as magnetic flux-freezing applies. A similar condition 
M/$ < (M/$) C rit = (2iiG 1,2 y 1 is required for uncondi- 
tional stability of an infinite uniform (in x,y) layer that is 
flatte ned along the z-direction o f a background magnetic 
field (jNakano fc Nakamural 1 19781 ). The various numerical 
values of (M/<&) cr it differ by small factors of orde r unity , 
and we adopt the result of iNakano fc Nakamural (|l978l ) 
since their model closely resembles the initial state in our 
calculation. 

Magnetic field strength measurements through the 
Zeeman effect reveal that the mass-to-flux r atios are 
clust e red a b out the critical value for collapse (|Crutcherl 
Il999t 12004 IShu et all Il999l ) and that there is also an 
approximate equipartition between the absolute values 
of gravitational e nergy and nonthermal (h e reafter , "tur - 
bulent" ) energy ([Myers fc Goodman! Il988l : iBasul l2000h . 
Measurements of polarized emission from dust grains, 
which reveal the field morphology, generally indicate 
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that the field in cloud cores is well ordered and not 
dominated by turbulent motions, with application of 
the Chandrasekhar-Fermi metho d yielding; mass -t o-flux 
ratios near the critical value dLai et al.l 200 ll. 20021; 
Crutcher. Nutter, fc Ward- Thompson! 2004; ICurran et all 



20041; iKirk. Ward-Thompson, fc Crutcherll2006l ) 

iMestel fc Spitzer] l|l956h pointed out that even if clouds 
are magnetically supported, ambipolar diffusion (resulting 
from ion-neutral slip) will cause the support to be lost and 
stars to form. More specifically, a medium with a subcrit- 
ical mass-to-flux ratio will still undergo a gravitationally- 
driven instability, occurring on the ambipolar diffusion time 
scale ra t her t h an the dynamical tim e scale (jLangerl 1 19781 : 
IZweibell 1 19981 ; ICiolek fe Basu] I2006T ). The length scale of 
the instability is essentially the Jeans scale in the limit 
of highly subcritical clouds (the same length scale as for 
highly supercritical fragmentation) but can be much larger 
when the mass-to-flu x ratio is close to the critical value 
l)Ciolek fc B asu 200£jj). Most nonlinear calculations of am- 
bipolar diffusion driven evolution have focused on a sin- 
gle axisymmetric core, but newer models focus on a frag- 
mentation process that results in the formation of multiple 
cores and somewhat irr e gular density and velocity structure. 
Ilndebetouw fc Zweibe ll (|2000h carried out a two-dimensional 
simulation of an infinitesimally-thin sheet threaded by an 
initially perpendicular magnetic field. Starting with slightly 
subcritical initial conditions, they followed the initial growth 
of mildly elongated fragments which occurred on a time scale 
intermediate between the dynamical time associated with 
supercritical collapse and the ambipolar diff usion time-scale 
associ ated with highly subcritical clouds. iBasu fc Ciolekl 
(2004) carried out two-dimensional simulations of a mag- 
netized sheet in the thin-disc approximation, which incor- 
porates a finite disc half-thickness Z consistent with hydro- 
static equilibrium and thereby includes the effect of mag- 
netic pressure. They studied a model which had an initially 
critical mass-to-flux ratio and another which was supercrit- 
ical by a factor of two. One of their main results was that 
the critical model had subsonic (maximum speed ~ 0.5c s , 
where c s is the isothermal sound speed) infall, while the 
decidedly supercritical cloud had infall speeds > lc s on 
scales ~ 0.1 pc from the core centres. This is a signifi- 
cant observationally-testable difference between dynamical 
(supercritical) fragmentation and ambipolar-diffusion reg- 
ulated (critical or subcritical) fragmentation. Yet another 
mode of fragmentation is the so-called turbulent fragmen- 
tation, which in fact corresp onds to collapse driven by a 
strong external comp ression. iLi fc Nakamural (|2004T ) and 
iNakamura fc Lil l|2005l ) have studied this process for a mag- 
netized sheet using the thin-disc approximation, and includ- 
ing the effect of ion-neutral friction. They find that a mildly 
subcritical cloud can undergo locally rapid ambipolar diffu- 
sion and form multiple fragments because of an initial large- 
scale highly supersonic compression wave. The core forma- 
tion occurs on a crossing time of the simulation box, w hich 
is related to the dynamical time. ILi fc Nakamural (|2004f ) es- 
timate that such a process can simultaneously maintain a 
relatively low efficiency of star fo rmation, as is required by 
observations (|Lada fc Ladall2003h . 

In this paper, we stu dy the three-dimensiona l exten sion 
of models such as those o f llndebetouw fc Zweibell (|2000l ) and 
IBasu fc Ciolekl (|2004T l. The self-consistent calculation of the 



vertical structure of the cloud allows us to test the predic- 
tions of two-dimensional models as well as to make some 
new predictions. We model clouds that are either decidedly 
supercritical or subcritical and study the evolution after the 
introduction of small-amplitude perturbations. The case of 
compression-induced collapse will be studied in a separate 
paper. We note that our three-dimensional model is not a 
cubic region but rather a flattened three-dimensional layer 
that is consistent with the expected settling of gas along 
the direction of the magnetic field. In reality, we believe 
that our modeled region would represent the dense midplane 
of a larger more turbulen t cloud. As demonstrated b y the 
one-dimensional models of iKudoh fc Ba su (2003, 2006|), the 
turbulent motions in a stratified magnetized cloud develop 
the largest amplitude motions in the outer low-density enve- 
lope, while maintaining transonic or subsonic motions near 
the midplane. We believe that the fragmentation process as 
modeled in this paper may proceed while long-lived turbu- 
lent motions continue on larger scales. 



2 NUMERICAL MODEL 
2.1 Basic Equations 

We solve the three-dimensional magnetohydrodynamic 
(MHD) equations including self-gravity and ambipolar diffu- 
sion, assuming that neutrals are much more numerous than 
ions: 



dp 

Tt +V 



Vp = -pV 



dv , 1„ 1 . „ _ , 

— + (v ■ X7)v = — Vp H j x B - VV>, 

at p cp 



OB 
~dt 



V x (v x B) + V x 



— (j xB)xB 
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V 2 ip = 4ttG P , 



P = c s p, 



(1) 

m 

(3) 

(4) 

(5) 
(6) 



where p is the density of neutral gas, p is the pressure, v is 
the velocity, B is the magnetic field, j is the electric cur- 
rent density, ip is the self-gravitating potential, and c s is the 
sound speed. Instead of solving a detailed energy equation, 
we assume isotherma l ity fo r each Lagrangian fluid particle 
j|Kudoh fc Basull2003l . l2006h : 



~dt 



dt 



— r— = —^r— + V ■ Vc s = 0. 



(7) 



For the neutral-ion collisio n time in equation ([3| and a sso- 
ciated quantities, we follow H asu fc Mouschoviaa (I1994T K so 
that 



= 1.4 



rrii + m„ 
pi{cw)in ' 



(8) 



where pi is the density of ions and {aw)i„ is the average 
collisional rate between ions of mass rrii and neutrals of mass 
m n . Here, we use typical values of HCO + -H2 collisions, for 
which {aw) in = 1.69 x 10~ 9 cm~ 3 s _1 and mi/m n = 14.4. 
We also assume that t he ion density Pi is determined by the 
approximate relation (|Elmegreenlll979l ; iNakanoll 19791 ) 
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pi = rriiK 



p/rrin 
10 5 cnT 



(9) 



where we assume K = 3x 10 3 cm 3 and k = 0.5 throughout 
this paper. 



2.2 Initial Conditions 

As an initial condition, we assume hydrostatic equilibrium 
of a self-gravitating one-dim ensional cloud along z-direction 
l|Kudoh fc Basij|2003ll2006t ). The hydrostatic equilibrium is 
calculated from equations 



dz 



= P9z 



%L = -4,Gp, 

dz 



P = CsP, 

subject to the boundary conditions 
fl,(* = 0) = 0, p(z = 0)=p o , p(z = 0) 



: Poc 2 s0 



(10) 

(11) 
(12) 

(13) 



where po and c„o are the initial density and sound speed 
at z = 0. If the initial sound speed (temperature) is uni- 
form throughout the region, we hav e the following analytic 
solution ps found bv lSpitzen (|l942h : 



ps(z) = posech (z/Hq) 
where 



Ho 



V^Gpo 



(14) 



(15) 



is the scale height. However, an isothermal molecular cloud 
is usually surrounded by warm material, such as neutral 
hydrogen gas. Hence, we assume the initial sound speed dis- 
tribution to be 



2.1/2 
c s0 T 2^ Cac 



2 \ 
CsO) 



1 + tanh 



(16) 



where we take c? ac — lOc^o, z c = 2Ho, and Zd = Q.lHo 
throughout the paper. By using this sound speed distribu- 
tion, we can solve equations (|10p -(ll2 p numerically. The ini- 
tial density distribution of the numerical solution shows that 
it is almost the same as Spitzer's solution for ^ z ^ z c . 

We also assume that the initial magnetic field is uniform 
along the z-direction: 



B z = Bq, B x = By 



0, 



(17) 



where Bq is constant. 

In this equilibri um sheet-like gas, we input a random ve- 
locity perturbation ( Mi vama. Narita. fc Havashilll987bh at 
each grid point: 



0.1c 3 oR 



O.lCsO-Rm, Dz 



0.0 



(18) 



where R m is a random number chosen uniformly from the 
range [-1,1]. The Rm'a for each of v x and v y are independent 
realizations. However, each model presented in this paper 
uses the same pair of realizations of R m for generating the 
initial perturbations. 



2.3 Numerical Parameters 

A set of fundamental units for this problem are c a o, Ho, and 
po- These yield a time unit to = Hq/c s o- The initial magnetic 
field strength introduces one dimensionless free parameter, 



Po 



8np 8np c a0 



a o 



IX 1 

°o 



(19) 



the ratio of gas to magnetic pressure at z = 0. 

In the sheet-like equilibrium cloud with a vertical mag- 
netic field, Po is related to the mass-to-flux ratio for Spitzer's 
self-gravitating cloud. The mass-to-flux ratio normalized to 
the critical value is 

Bo 

where 



PS 



(20) 



Ss = 



psdz = 2p H 



(21) 



is the column density of Spitzer's self-gravitating cloud. 
Therefore, 



Po 



2 

PS- 



(22) 



Although the initial cloud we used is not exactly the same 
as the Spitzer cloud, Po is a good indicator to whether or 
not the magnetic field can prevent gravitational instability 
ijNakano fc Nakamuralll978r i. 

Dimensional values of all quantities can be found 
through a choice of po and c s o- For example, for c s o = 0.2 
km s _1 and no = po/m n = 10 4 cm -3 , we get Ho = 0.05 pc, 
t = 2.5 x 10 5 yr, and B = 20 p,G if Po = 1. 



2.4 Numerical Technique 

In order to so lve the equat i ons numerically, we use the 
CIP method (|Yabe fc Aokil Il99ll : lYabe. Xiao, fc Utsumil 
l200lT ) for equations flTJl, © and ©, and the method 
of characterist ics -cons trained transport (MOCCT; 
IStone fc Normanl l| 19921 ) ) for equation (O, including 
an explicit integration of the ambipolar diffusion term. 
The combination of the CIP and MOCCT me th ods is 
Kudoh. Matsumoto. fc Shibatal (|l999h 



sum marized 

and lOgataet al.l d2004h . It includes the CCUP method 
l|Yabe fc Wand 1 199 £ for the calculation of gas pressure, 



order to get more numerically stable r esults. The numeri cal 
code in this paper is based on that of lOeata et aD l|2004 ). 

In this paper, the ambipolar diffusion term is only in- 
cluded when the density is greater than a certain value, 
per- We let p cr = 0.3po both for numerical convenience 
and due to the physical idea that the low density region 
is affected by external ultraviolet radiation so that the 
ionization fraction becomes large, i.e. r n i becomes small 
IjCiolek fc Mouschoviasl Il995i ). Under this assumption, the 
upper atmosphere of the sheet-like cloud is not affected by 
ambipolar diffusion. This simple assumption helps to avoid 
very small time-steps due to the low density region in order 
to maintain stability of the explicit numerical scheme. 

We used a mirror-symmetric boundary condition at 
z = and periodic boundaries in the x and y-directions. 
At the upper boundary at z = z out = 4Ho, we also used 
a mirror-symmetric boundary except when we solve the 
gravitational potential. This symmetric condition is just 
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Figure 1. The time evolution of the density at the location 
where the density reaches its maximum value in each simulation. 
Each line shows a case of different /3q- 

for numerical convenience. However, because the results 
we show later in this paper are consistent with previous 
two-dimensional simulations, we believe that the boundary 
conditions do not affect the result significantly. The Pois- 
son equation <(Sj is solved by the Greens function method 
to compute the gravitational kernels in z-direction, along 
with a Fourier transform method in the x and ^-directions 
|Mivama. Narita. fc Havashilll987bl ). This method of solv- 
ing the Poisson equation allows us to find the gravitational 
potential of a vertically isolated cloud within \z\ < z out . 

The computational region is \x\, \y\ ^ 8nHo and ^ 
z ^ 4Hq. The number of grid points for each direc- 
tion is (N x ,N y ,N z ) = (64,64,40). Since the most un- 
stable wavelength for no magnet ic field is about 4irHo 
l|Mivama. Narita. fc Havashilll987al ). we have 16 grid points 
within this wavelength. We have also 10 grid points within 
the scale height of the initial cloud in the z-direction. While 
this is not a high-resolution simulation, we believe that we 
have the minimum number of grid points to study the grav- 
itati onal instability, es pecially by using the code based on 
CIP (|Ogata et alj2004h . The maximum computational time, 
which occurs for the case of the subcritical cloud, is about 85 
hours of CPU time using a single processor of the VPP5000 
in the National Astronomical Observatory of Japan. 



3 RESULTS 

Figure 1 shows the time evolution of the density at the lo- 
cation where the density reaches its maximum value in each 
simulation. The simulations are stopped when the maximum 
density is about 30po- Each line shows a case of different f3o- 
When /So is 100 or 4, the magnetic field is not strong enough 
to suppress the self-gravitational instability of the cloud. 
In these cases, the density evolves rapidly, over the sound- 
crossing time of the most unstable wavelength (~ 4ivHo)- 
However, when /3o = 0.25, the cloud is self-gravitationally 
stable unless ion-neutral slip is present. Therefore, the den- 
sity evolves gradually over the diffusion time of the magnetic 
field. According to the two-dimensional linear analysis by 
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Figure 2. The logarithmic density contours for /3q = 0.25 at t = 
150. Arrows show velocity vectors on each cross section. Upper 
panel shows the cross section at z = 0, and the lower panel shows 
the cross section at y = 4.3. 

ICiolek fc Basul (|2006h . the evolutionary time scale of a sig- 
nificantly subcritical cloud is about ten times longer than the 
dynamical time, for a standard ionization fraction, as used 
here. Our numerical result is consistent with their analysis. 

Figure 2, Figure 3, and Figure 4 show the logarithmic 
density contours for f3o — 0.25 at t = 150, /3o = 4 at t — 15.3, 
and f3o ~ 100 at t = 11.1 respectively. Each upper panel 
shows the cross section at z — 0, and the lower panel shows 
the cross section at y — 4.3 for /3q = 0.25, y = 5.1 for /3o = 4, 
and y = 5.1 for /3o = 100 respectively. The values of y for 
the lower panels are chosen so that the vertical cut passes 
through at least one dense core. (In these numerical simu- 
lations, we use the term "core" to refer to the region where 
the density is greater than the mean background density by 
about a factor of 3.) The size of cores for /3o = 4 is bigger 
than that of po = 100. The size becomes smaller again when 
the magnetic field is stronger than critical (/3q = 0.25). This 
result is consistent with th e two-dimensional linear analy- 
sis of ICiolek fc Basul (2006), who found a hybrid mode for 
critical or mildly supercritical clouds in which the combined 
effect of field- line dragging and magnetic restoring forces en- 
force a larger than usual fragmentation scale. Arrows show 
velocity vectors on each cross section. Maximum velocities 
become supersonic for /3o = 4 and /3o = 100, but remain 
subsonic for 0o = 0.25. This is also c onsistent with the two - 
dimensional numerical simulations of lBasu fc Ciolekl (|2004h . 

Figure 5 and Figure 6 shows the close up views of the 
density contours around cores for flo = 0.25 and f3o — 4 
respectively. Magnetic field lines near cores are also plotted 
in three-dimensional space. When /3q = 0.25, the neutral gas 
has to slip through the field lines to make a gravitationally 
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Figure 3. The logarithmic density contours for /3q = 4 at t = 
15.3. Arrows show velocity vectors on each cross section. Upper 
panel shows the cross section at z = 0, and the lower panel shows 
the cross section at y = 5.1. 




Figure 4. The logarithmic density contours for (3q = 100 at 
t = 11.1. Arrows show velocity vectors on each cross section. 
Upper panel shows the cross section at z = 0, and the lower 
panel shows the cross section at y = 5.1. 
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Figure 5. The close up views of the density contours around 
cores for /3q = 0.25. Magnetic field lines near cores are also plotted 
in three-dimensional space. 




Figure 6. The close up views of the density contours around 
cores for /3o = 4. Magnetic field lines near cores are also plotted 
in three-dimensional space. 

bound core. Therefore, the field lines are not so deformed in 
the case of fio = 0.25, in contrast to those of /3o = 4. 

Figure 7 and Figure 8 show the logarithmic contours of 
plasma j3 for /3q = 0.25 at t = 150 and /3o — 4 at t = 15.3, re- 
spectively. When Po = 0.25, the plasma (5 in cores is greater 
than in the surroundings. This is because the mass-to-flux 
ratio (and therefore (S) has to increase in order for the core to 
become gravitationally unstable. The maximum f3 is larger 
than 1 at the centre of a core, which means that the centre of 
the core is approximately supercritical. In contrast to this, 
the plasma (3 in cores is slightly lower than the surround- 
ings when Po = 4. In this case, the magnetic field is swept 
up by the contracting cloud before ion-neutral slip works 
efficiently. If hydrostatic equilibrium along the z-direction 
is exactly satisfied, the plasma f3 would remain constant in 
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Figure 9. The densities, plasma /3, and v x along rr-axes for lines 
that cut through the cores shown in Figure 5 and Figure 6. The 
left panel shows the core for /3o = 0.25. The right panel shows the 
core for j3o = 4. 



Figure 7. The logarithmic contours of plasma /3 for (3q = 0.25 
at t = 150. Upper panel shows the cross section at z = 0, and the 
lower panel shows the cross section at y = 4.3. 
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Figure 8. The logarithmic contours of plasma f3 for /3rj = 4 at 
t = 15.3. Upper panel shows the cross section at z = 0, and the 
lower panel shows the cross section at y = 5.1. 
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Figure 10. The densities, plasma (3, and v z along 2-axes for 
lines that cut through the cores shown in Figure 5 and Figure 
6. The left panel shows the core for /3rj = 0.25. The right panel 
shows the core for /3q = 4. 



time and space. The slightly lower values of (5 in cores are 
probably caused by the modestly nonequilibrium state along 
2 during the evolution. 

Figure 9 shows the densities, plasma (3, and v x along 
:r-axes for lines that cut through the cores shown in Figure 
5 and Figure 6. The left panel shows the core for /3q = 0.25. 
The right panel shows the core for (3q = 4. Figure 9 shows 
that the plasma (3 in the core is higher than the surroundings 
for Po = 0.25 and lower for /3o = 4. The infall velocity for 
Po = 4 shows supersonic values, while the velocity for /3q = 
0.25 is subsonic. These vel ocities are comparable to those in 
Figure 2 and Figure 3 of lBasu fc Ciolekl (|2004 ). Figure 10 
shows the densities, plasma /3, and v z along z-axes for lines 
that cut through the same cores. It also shows that the infall 
velocity for f3o = 4 reaches mildly supersonic values, while 
the velocity for f3o = 0.25 is subsonic. 

Figure 11 shows the relation between density and mag- 
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Figure 11. Open circles show the magnetic field strength as a 
function of density along z = at t = 150 for the model with 
fa = 0.25. The strength of the magnetic field is normalized by 
yj 8npQC^ . Filled circles are the same for /3o = 4, at t = 15.3. The 
blue line shows the evolutionary track of the point at which the 
density achieves its maximum value for the model with /3o = 0.25. 
The red line is the same for /3q = 4. 



netic field on the plane z — 0. T he strength of the magnetic 
field is normalized by -\/87rpo c 2 sQ . Open circles show the mag- 
netic field strength as a function of density along z = at 
t = 150 for the model with fa = 0.25 (see Fig. 2). Filled 
circles are the same for fa = 4, at t = 15.3 (see Fig. 3). The 
blue line shows the evolutionary track of the point at which 
the density achieves its maximum value for the model with 
fa — 0.25. The red line is the same but for fa = 4. This fig- 
ure shows that the snapshot of the relation between density 
and magnetic field at different spatial points in the midplane 
of the cloud overlaps with the evolutionary track of an in- 
dividual core . The dashed line shows B oc p ' 5 . When the 
density becomes large, each relation approximately tends to 
B oc p 0,5 . In the case of fa — 0.25, the relation shows that 
core initially evolves to greater density without increasing 
the magnetic field strength. This is caused by the slip of 
neutral gas through the magnetic field during the subcriti- 
cal phase of evolution. 



4 CONCLUSIONS AND DISCUSSION 

We have studied fragmentation of a sheet-like self- 
gravitating cloud by three-dimensional MHD simulations. 
The main results are as follows. 

• We confirmed that in the case of an initially subcritical 
cloud (fa — 0.25), cores developed gradually over an am- 
bipolar diffusion time, while the cores in an initially super- 
critical cloud (fa — 4 or fa — 100) developed in a dynamical 
time. 

• The infall speed on to cores is subsonic in the case of 
an initially subcritical cloud, while there is extended super- 
sonic infall in the case of an initially supercritical cloud. This 
is consi stent with the result of the two-dimensional simula- 
tions of lBasu fc Ciolekl(|2004h . In our three-dimensional sim- 



ulations, we also find that the ^-component of the velocity 
follows the same pattern. 

• The size of cores for mildly supercritical cloud (fa = 4) 
is bigger than that of highly supercritical cloud (fa = 
100). The size becomes smaller again when the magnetic 
field is stronger than the critical (fa — 0.25). This re- 
su lt is consistent with the two-dimensional linear analysis 
of ICiolek fc Basul (|2006n . 

• When the cloud is initially subcritical (fa = 0.25), the 
plasma j3 in cores is greater than in the surroundings. In 
contrast to this, the plasma /3 in cores is slightly lower than 
the surroundings when the cloud is initially supercritical 
(fa = 4). The latter result is probably caused by the mod- 
estly nonequilibrium state along z during the evolution. 

• In the B — p plane, the snapshot of the relation between 
magnetic field strength (B) and density (p) at different spa- 
tial points of the cloud overlaps with the evolutionary track 
of an individual core. When the density becomes large, each 
relation approximately tends to B oc p ' 5 . 

Our simulation is the first fully three-dimensional sim- 
ulation to study the role of magnetic fields and ion-neutral 
friction in fragmentation. In this paper, we concentrated on 
the effect of initially small perturbations, partly as a way to 
compare with established predictions of linear theory. Our 
models also serve as a guide to understand fragmentation oc- 
curring exclusively in dense subregions of clouds that contain 
only subsonic or transonic motions. Real molecular clouds 
certainly contain supersonic turbulence, at least in their low- 
density envelopes, as is observed through large line-widths 
of emission lines from relatively low-density tracers. The ad- 
ditional effect of supersonic turbulence on three-dimensional 
fragmentation with magnetic fields and ion-neutral friction 
will be studied in an upcoming paper. 
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